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ABSTRACT 

Detectable flux is visible in the Lya and Ly/3 troughs of the highest redshift (z — 6.42) quasar found 
to date, SDSS J1148+5251. This has previously been interpreted as continuum contamination from 
an interloper galaxy a,t z — 4.94. We examine the Ly7 trough of SDSS J1148+5251 and show that this 
interpretation is untenable: the spectrum does not show the continuum break in a z = 4.94 galaxy 
expected from absorption by the intervening Lya forest. Therefore, flux must be leaking through 
at least one of the troughs from the quasar itself. Contrary to previous claims, the flux ratios in 
the Lya and Ly/3 troughs are consistent with pure transmission. From the Ly7 trough, we place an 
upper bound on the effective Lya optical depth at z ~ 6.2 of Tee < 14.3 (2a). This implies a highly 
ionized IGM along this line of sight and significant cosmic variance in the transition toward complete 
Gunn-Peterson absorption. Detailed study of the observed transmission features will shed light on 
this era. 

Subject headings: cosmology: theory - early universe - galaxies: formation -high-redshift — galaxies: 
quasars: general-absorption lines 



1. INTRODUCTION 

Despite spectacular progress in recent years, the reion- 
ization history of the universe remains shrouded in mys- 
tery. The high electron scattering optical depth observed 
by WMAP indicates that reionization may have begun 
as early as z ~ 20 (Kogut et al 2003; Spergel et al 2003). 
On the other hand, spectra of high redshift quasars indi- 
cate a sharp change in the neutral hydrogen fraction and 
ionizing background at z ~ 6, suggesting this era marks 
the end of the reionization epoch (Becker et al 2001, Fan 
et al 2002). Taken in tandem, these observations suggest 
a highly complex reionization history. 

A crucial missing piece of the puzzle is the precise ion- 
ization state of the IGM at z = 6. The large optical 
depth of the IGM to hydrogen Lyman transitions im- 
plies that Gunn & Peterson (1965) absorption can (at 
best) constrain the volume and mass-weighted neutral 
fractions to be xhi,v > 10~^, xm,M > 10~^ respectively 
(Fan et al 2002). 'it is plausible that the IGM is stiU 
highly ionized at z = 6, and with full reionization occur- 
ring much earlier. On the other hand, modeling of the 
spectral regions around two quasars at z > 6.2 may indi- 
cate larger neutral fractions (xhi ^0.2; Wyithe & Loeb 
2004a; Mesinger & Haiman 2004). At stake is the very 
nature of reionization: if these claims are correct, then 
reionization must be an extremely rapid process, since 
the IGM is known to be highly ionized (xhi < 10~^) by 
z = 5.9 along all observed lines of sight (Fan et al 2003). 

It is thus worth considering the IGM absorption in 
more detail. In particular, a completely dark Gunn- 
Peterson trough at z ~ 6 is not necessarily universal. 
Both the Lya and Ly/3 troughs of the highest redshift 
quasar found to date, SDSS J1148-f5251, contain de- 
tectable flux (White et al 2003). In this Letter, we show 
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that its previous interpretation as continuum contamina- 
tion from an interloper galaxy at z — 4.94 is inconsistent 
with the observed flux ratios blueward and redward of 
the expected continuum break due to Lya forest absorp- 
tion. Therefore, some of the flux is true transmission 
due to holes in the high redshift Lya forest. We derive 
our strongest contraints from the Ly7 trough, which was 
ignored in previous analyses and should be relatively un- 
contaminated by flux from an interloper. The presence 
of flux transmission implies that either the IGM is still 
highly ionized at z = 6 or that there is significant cosmic 
variance in the reionization epoch along different lines of 
sight. 

2. THE LYMAN GAMMA TROUGH IN SDSS J1148-I-5251 

Detectable flux and a network of transmission fea- 
tures is seen in both the Lya and Ly/3 troughs of SDSS 
J1148-I-5251. White et al (2003) interpreted this to 
be continuum contamination from interloper galaxies at 
z = 4.94 for two reasons: (i) strong C IV absorption 
features appear at z = 4.9, so the apparent strong Ly/3 
spikes at z — 6.03, 6.06 could just be Lya emission from 
z w 4.94. (ii) The flux ratios of the troughs appear incon- 
sistent: there is too little light in the Ly/3 trough given 
the Lya transmission. Although the a priori probability 
of such an interloper is small, it increases if the interloper 
also lenses SDSS J1148-f 5251, boosting its probability of 
detection. 

By examining the Ly7 trough of SDSS J1148+5251, 
we argue that such an interpretation is untenable. We 
first note that if residual flux in any of the Lya,/3,7 
troughs can be attributed to the quasar rather than an 
interloper galaxy, the IGM must still be highly ionized 
along the line of sight. Consider the boundary of the 
quasar proximity zone at zhii ~ 6.33, where the Lya, 
Ly/3 absorption increase rapidly. There is also a sharp 
jump in absorption at the boundary of the Ly^ trough, 
strongly suggesting that the transmission features redward 
of the Ly^ trough are genuine. Even if we exclude the 
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sharp transmission spike at A = 7205 A (putatively a Lya 
emission hne from a galaxy at z = 4.94), the residual flux 
changes from F_2o{z-y = 6.33-6.40) = 29.0±1.5 (within 
the quasar proximity zone) to F2o{zj = 6.25 — 6.32) = 
4.3±1.4 (outside the quasar proximity zone), where -F_2o 
is the flux F in units of 10~^''ergs~^ cm~^ A~^. It would 
be remarkable if the protocluster at z w 5 lined up so as 
to exactly coincide with the quasar proximity zone in 
Ly7: it is much more likely that the sudden change is 
due to genuine Ly7 absorption. This suggests the flux 
seen in the Ly/3 forest at z > 5.95, before the onset of 
Ly-f absorption, represents true transmission. 

It is also crucial that there is detectable flux in the 
Ly7 trough. Consider the wavelength stretch corre- 
sponding to z = 6.17 — 6.32 in the various absorption 
troughs. This lies outside the quasar's apparent re- 
gion of influence and ends where the Ly7 trough be- 
comes contaminated by LyJ and higher order absorption 
(A = 6962 A, corresponding to z-y = 6.16). White et 
al (2003) suggest that the Lya, Ly/3 troughs are con- 
taminated by continuum emission from the interloper 
galaxy so that detected flux is not significant. To be con- 
servative, let us ignore the transmitted flux in the Ly/3 
trough, which is potentially also contaminated by emis- 
sion line features from the z « 5 proto-cluster. We find 
that F_2o(Lya)= 3.0 ± 0.9 while i^_2o(Ly7)= 4.9 ± 0.9; 
these regions correspond to A = 1467 — 1498 A and 
A = 1174 — 1199 A in the rest frame of an interloper 
galaxy at z = 4.94. These flux ratios are strongly incon- 
sistent with the expected continuum of a z k, b galaxy, 
which should have a strong break due to the intervening 
Lya forest. Songaila & Cowie (2002) obtain a flux sup- 
pression factor Tq, = 0.14 ±0.03 for the Lya forest in this 
redshift and wavelength interval (Becker et al 2003 find 
Tq = 0.11). Even if the flux in the Lya trough is entirely 
due to continuum contamination from the interloper, the 
implied continuum contribution to the Ly7 trough would 
be F_2o = 0.4(T„/0.14) < i^_2o(Ly7), which lies weU 
within the noise. Thus, almost all the observed flux in the 
Ly7 trough must be genuine transmitted flux from the 
quasar. Because it is almost completely free from contin- 
uum contamination, the Ly7 trough places the strongest 
constraint on the IGM ionization state along this line of 
sight. 

Because there are indisputably C IV absorbers at 
z = 4.9, the putative protocluster may conceivably have 
another effect: it could ionize the IGM at z = 5, cre- 
ating a transparent Lya window there. However, the 
proximity zone of the interloper galaxy itself is too small. 
We need to highly ionize the IGM on comoving length- 
scales L w 55(Az/0.1) Mpc at z = 5. If the flux 
in the Lya trough is interpreted as the continuum of 
the interloper, then F_2o(A = 1450A) = 3.9. If we 
take F{\ = 1450l)/i^(A 900l) = 4.6 as is ob- 
served in z = 3.4 LBGs (Steidel et al 2001) for the 
strongest ionizing continuum possible (corresponding to 
a ~ 100% escape fraction for ionizing photons), then we 
infer L^lVRyd) — 1.4 x lO^^erg s~^ cm~^ Hz~^, roughly 
implying a star formation rate SFR ~ 14MQyr~^ for 
a Salpeter IMF. Assuming the mean background ioniz- 
ing rate in units of 10~^^s~^ to be r_i2(z = 5) w 0.15 
(e.g., Fig 2 of Fan et al 2002), the local radiation field be- 
comes comparable to the metagalactic ionizing at r « 1.3 



Mpc comoving, far too small to be of interest. Lyman- 
break galaxies (LBGs) at z « 3 — 4 also show excess 
Lya transmission on lengthscales i w 0.5 Mpc comoving 
(Adelberger et al 2003), roughly corresponding to the 
distance a 600kms~^ wind would travel on the ~ 300 
Myr star formation timescale of LBGs. It is unlikely a 
z = 5 galaxy would affect a significantly larger volume. 
Galaxies near (but not along) the line of sight could con- 
tribute additional ionizing photons. Note, however, that 
we require a factor ~ 5 increase in Tq, over a lengthscale 
~ 10 times larger than the comoving correlation length 
for highly biased galaxies. There is no evidence for such 
strong fluctuations in transmission at this redshift: typ- 
ically (JTa/Ta ~ 0.3 (Songaila & Cowie 2002). In any 
case, an anomalously large protocluster should be easily 
visible in narrowband optical searches. 

3. IMPLICATIONS FOR THE IGM AT Z = 6 

Having established the reality of residual flux from 
SDSS J1148-I-5251, we will now consider its implications 
for reionization at z ~ 6. We must first estimate the 
effective optical depths Tcff in each transition. We begin 
by summarizing our methodology and the error budget 
(which is frequently underestimated). 

The errors for Ly7 transmission T-y = Ttot/(7Q,T^) 
must include not only pixel noise but also the cosmic vari- 
ance in foreground transmission in Lya and Ly/3. Pre- 
vious analyses have often neglected this additional scat- 
ter. It has been measured in quasar samples; we use the 
values in Table 2 of Songaila & Cowie (2002), who tabu- 
late the measured {Ta),a'Ta for 6 redshift bins between 
z = 4.1 — 5.5. However, these use fixed wavelength in- 
tervals, and (TTa obviously depends on the corresponding 
comoving length L. We cannot simply assume Poisson 
statistics, because long wavelength modes could domi- 
nate the variance. To estimate (Ttq we follow the ansatz 
of Lidz et al (2002), who argue that the transmission 
power spectrum takes the shape (though not the nor- 
malization) of the linear mass power spectrum on large 
scales, yielding: 
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where (T„ is the noise per pixel, N is the number of pixels, 
and Pf{k) — Bexp{—ak^) {dk/2TT)kP^s,ss(k), where 

a K kj kj is the Jeans wavenumber, and B is nor- 
malized to the observed CTa for some observed stretch 
Lobs- The first term in equation typically dominates 
by an order of magnitude. For the large lengthscales 
of interest, Pf{k) is fairly flat, and we have approxi- 
mately (T^^ oc L~^, as expected for a white noise power 
spectrum. Note that this ansatz for the transmission 
power spectrum Pf(k) assumes a homogeneous ionizing 
background F; if F exhibits small-scale fluctuations, cttq 
could be larger. However, there is no cosmic variance 
in Tfj, because we have observed the corresponding Lya 
transmission (at Zq, = 5.9) along the same line of sight. 
The main uncertainty comes from the optical depth con- 
version between different lines (see below) . 

Another crucial item is uncertainty in the quasar con- 
tinuum. This is particularly important in comparing rel- 
ative absorption between the different troughs. The Lya 
line is much broader than the Ly/3 and Ly7 lines and gen- 
erally spills over into the observed trough. Uncertainty 
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TABLE 1 

Absorption Optical Depths from z = 6.17 - 6.32 in SDSS ,11148+5251 



Trough 


Flux'' 


Continuum^"'^ 


Total Optical Depth 


Line Optical Depth'^ 


rc(/3 = 2)d 


T<,(/3 = 3)d 


Lya 
Ly/3 
Ly7 


3.0 ±0.9 
9.0 ± 1.1 
4.9 ±0.9 


(2560,1730,1690) 
(2140,1870,1710) 
(2070,1910,1710) 


(6.8,6.4,6.4) ±0.3 
(5.5,5.3,5.2) ±0.1 
(6.1,6.0,5.9) ±0.3 


(6.8,6.4,6.4) ±0.3 
(3.1,3.0,2.9) ±0.3 
(2.5,2.4,2.3) ±0.3 


(6.8,6.4,6.4) ±0.3 
(8.2,7.8,7.6) ±0.9 
(11.5, 11.2, 10.7) ± 1.4 


(6.8,6.4,6.4) ±0.3 
(6.7,6.4,6.2) ±0.7 
(8.4,8.2,7.8) ± 1.0 



In units of lO^^^crgs"^ cm"^A"\ 

^Derived from the LBQS composite spectrum, and the Telfer et al 2002 spectral indices for radio-quiet (qeuv — —1.57) and radio-loud (qeuv — 
— 1.96) quasars respectively. Subsequent triplets for the calculated optical depths reflect this ordering. 

^ After subtracting off the contribution from foreground absorption. For Ly/3 trough; Tfg — Ta{z — 5.1) — 2.38 ± 0.32; for Ly7 trough: Tfg — 
T^{z = 4.8) + Tp(z = 5.9) = (1.95 ± 0.25) + (1.59 ± 0.02) = 3.54 ± 0.25. 

*^The effective equivalent Lya optical depth inferred from the nonlinear conversion described in the text, assuming t(A) (X A^. 



in the strength of the hne creates significant uncertainty 
in the inferred Tq. A simple way to quantify this is to 
compare the optical depths inferred assuming the com- 
posite spectrum from the Large Bright Quasar Survey 
(LBQS; Brotherton et al 2001) and the far-UV quasar 
power law spectrum of Telfer et al (2002). The former 
includes all the emission line structure, while the latter 
assumes a pure power law; the difference between the two 
therefore captures the uncertainty in emission line con- 
tribution. We normalize the continuum to the observed 
flux at ~ 1290 A. 

Table 1 lists the measured effective line optical depths 
Ti = Teff(Lya) and their associated la errors. The results 
for the Lya and Ly/3 line optical depths are comparable 
to those of White et al. (2003). They argued that the 
pair are incompatible with each other, because oc Xifi 
(where fi is the oscillator strength of line i), implying 
Ta/r/s = 6.24 and Ta/r-y = 17.93. However, these rela- 
tions are only true at fixed density, and only apply if the 
IGM is homogeneous. For the large comoving lengths 
which we are considering, the flux transmission comes 
from a variety of densities: 

(T,) = (exp (-Teff,,)) = J exp [-r,(A)] P{A)dA, (2) 

where i^(A) is the probability distribution of overden- 
sities A = p/p, n oc X^f^{l + z)*-'^A^a{T)/T (e.g., Hui 
& Gnedin 1997) in the optically thin limit, and a{T) is 
the recombination coefficient. We use the form of P(A) 
given by Miralda-Escude, Haehnelt & Rees (2000), which 
is a good fit to numerical simulations. If we assume an 
equation of state T oc A"'' (where 7 ~ 0-1) and a fluc- 
tuating radiation field whose amplitude may be density 
dependent, T = FoA*, then n = A{z){l + z)^-^A'^ , where 
/3 = 2 — 0.77 — We can solve for the normalization 
constant A{z) by demanding {Ti{A)) — Tohs{z)- Note 
that A{z) oc f^X^/To. 

The integral in equation (|2Jl can be evaluated by the 
method of steepest descents (Songaila & Cowie 2002, 
Songaila 2004): 

reff = -AV(i+3/./4) ^ ^^^ln(A) + const. (3) 

The leading term implies that Toft oc A'^-^ ioi P — 2 (cor- 
responding to a uniform radiation field and an isother- 
mal equation of state). Thus, in an inhomogeneous uni- 
verse the optical depth increases more slowly than lin- 
early with the oscillator strength; conversely, it increases 



more slowly with a weaker radiation field. Evaluating 
equation Q numerically, we find that Tq/t^ « 3, and 
Ta/Tj w 5 — 6, with weak dependence on redshift, the 
equation of state, and F. This is because transmission is 
dominated by rare voids, and the primary effect of de- 
creasing fi is to increase the range of densities sampled 
by the line. 

This behaviour becomes even more marked if F is not 
uniform, which is certainly the case before reionization 
is complete. As a fiducial case, consider a uniform ra- 
diation field F = 0.05 at z = 6.15, which produces 
a Lya effective optical depth Tq, = 7. In this case, 
{Ta/Tf3,Ta/Tj) = (2.7,4.9). Now let us consider vari- 
ous scenarios that would produce the same Tq, but dif- 
ferent {Tp,T-y). Suppose the optical depth r oc A'' has 
P > 2, which could, for example, mimic self-shielding 
in overdense regions. In this case, from equation Q, 
the optical depth increases even more slowly with oscil- 
lator strength: for f] — 3, {Ta / Tfs , Ta / t^) = (2.1,3.4). 
Only if the optical depth is independent of overdensity 
(/? = 0) will we recover the linear scaling oc A. Fluc- 
tuations in the radiation field that are uncorrelated with 
density fluctuations (so that r oc A^) will produce simi- 
lar effects. For example, a radiation field with a lognor- 
mal probability distribution (F,crinr) — (0.02,1) yields 
(r„/r^,T„/rT,) = (1.9,2.9). 

Thus, the fluctuating density and radiation flelds in- 
troduce considerable uncertainty in the relations between 
TaiTjs, and Tj. The effective optical depth is often used 
to infer a;Hi,v a-nd xbi.m in the IGM, a procedure that 
we caution is fraught with uncertainty. In the case of 
a^HLM, it is almost meaningless. It is easy to see why: we 
can infer (xhi) oc (r) from (e^^) only in the limit where 
r ^ 1 and (e^"^) w 1 — (r); otherwise, our ignorance of 
the full probability distribution P{t) means that a wide 
variety of (xhi) would be consistent with a given observed 
(e^"^). As an illustration, we show in Figure 1 the loga- 
rithmic integrand j/AP(A) oc (y), for various quantities 
y; all the curves have been normalized to have unit area. 
Since (e""^) is heavily weighted toward voids and (a;Hi) is 
weighted toward overdense regions, estimating one from 
the other is extremely model-dependent. In particular, 
a^Hi.M is heavily weighted toward large overdensities, and 
the quasar spectra essentially leave it unconstrained. For 
instance, the neutral fraction could be considerably lower 
in overdense regions without affecting the transmitted 
flux. 
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Fig. 1. — The logarithmic integrand yAP{A) on (y) as a function 
of A for (Tq), (Tp), (T-y), (xhi,v)i (^^hi.m)- We assume a uniform 
ionizing background F = 0.04 at z = 6.15 and an isothermal equa- 
tion of state. The integrand for (xhi^m) peaks at high overdensities 
to the right of the plot. As the overlap integral between the dif- 
ferent Lyman series lines drops, the uncertainty in inferring one 
from the other increases. More importantly, there is considerable 
uncertainty in inferring (xhi) from quasar spectra; in particular, 
estimates for (xhi,m) are almost meaningless. 



With these caveats, we hst Tq, derived from each of the 
three transitions in the rightmost columns of Table 1 for 
(3 = 2, 3. The flux ratios in all 3 troughs are entirely 
consistent with flux transmission from the quasar with- 
out any contamination from an interloper: the equivalent 
Lya optical depths all lie within 1 — 2cr of one another. 
The apparent inconsistency in White et al (2003) ap- 
pears simply because they did not integrate over the IGM 
density distribution. Again, we caution that the optical 
depth conversion between different transitions has large 
uncertainties (not reflected in the error bars) due to un- 
certainty in the probability distribution of optical depths 
P(t). The true effective optical depth inferred from the 
Lya and Ly/3 troughs could also be somewhat higher 
if there is some continuum contribution from an inter- 
loper, which we cannot completely rule out. The most 
stringent optical depth constraint comes from the Ly7 
trough, which implies Tq < 14.3 (2cr), with a most likely 
value r„ w 6-10. By contrast, in SDSS J1030+0524, 
the ItT {2a) lower limit to the optical depth in the Ly/3 
trough is Tcff > 11.1 (9.9); the Ly-f trough yields similar 
constraints. 

4. DISCUSSION 

In this Letter^ we have argued that the residual flux in 
the Lya,/3, and 7 troughs of SDSS J1148+5251 is not 
due to an interloper galaxy but represents true transmis- 
sion in the z > 6 IGM. This places an upper bound on 
the effective Lya optical depth of Tcs < 14.3 (2(7), im- 
plying that the IGM is still highly ionized at z ~ 6.3. 
It has been argued that the size of the H II regions of 
the two highest redshift quasars i?Hii ~ 4.5 Mpc im- 
ply xm > 0.1 in this range (Wyithe & Loeb 2004a): 



i?Hii « 7a;jjy^(tagc/10'^ yr)^/^ Mpc, where tago is the hfe- 
time of the quasar. This is of course strongly dependent 
on the assumed template spectrum; if the escape frac- 
tion of ionizing photons in high-redshift quasars is some- 
how smaller, the constraint weakens. The quasars could 
also be lensed and hence intrinsically fainter (Haiman 
& Gen 2002), though follow-up HST observations have 
failed to detect multiple images in either quasar (White 
et al 2004 in preparation). For our purpose, we note that 
H II region radius is determined by the mass-weighted 
neutral fraction xhi.m- As we showed in this is 
very poorly constrained by Toff, because voids contain- 
ing only a small fraction of the mass dominate the trans- 
mission. For instance, if xhi.m ~ 0.1 of the baryons 
are in neutral self-shielded halos, we could still have 
xm.y ~ xhim/^ 0.1/200 ~ 5 x 10"'*, which could 
easily be accommodated by our results. 

A second argument in favor of xm > 0.1 is an indirect 
detection of the Gunn-Peterson damping wing (Mesinger 
& Haiman 2004). There is a stretch at the boundary of 
the H II region of SDSS J1030-I-0524 with Ly/3 transmis- 
sion but no Lya transmission, implying 6 < < 11. 
The absence of any Lya transmission from low-density 
voids in this segment implies a source of smooth opacity, 
attributed to a Lya damping wing that requires a large 
optical depth t > 10^. However, SDSS J1148-h5251 con- 
tains no such transition region. We also caution that 
the statistical significance of such transition regions is 
still unclear. For instance, in SDSS J1148-f5251 there 
is a stretch from z = 5.95 — 6.0 with Ly/3 transmis- 
sion {F^20 — 37.9 ± 1.7) but no significant Lya flux 
(F_20 — 2.5 ± 1.7). Again, the optical depth ratios are 
consistent with pure flux transmission (tq = 6.6 ± 0.7, 
Tfj = 2.0 ± 0.4 Tcff^Q « 5.5 ± 1.2). This stretch is 
~ 2.5 times longer than the Az = 0.02 zone seen in SDSS 
J1030-I-0524, over which even damping wing absorption 
would change substantially; in any case, we have argued 
that the IGM is highly ionized along this stretch. Such 
regions may therefore be fairly generic and not indicative 
of a damping wing. More detailed analysis, incorporat- 
ing a pixel-by-pixel analysis, variance in the foreground 
Lya forest, and better theoretical modeling of the fluc- 
tuating radiation and density fields in realistic models of 
reionization, would help shed light on this issue. 

Nonetheless, if the region around SDSS J1030-f0524 is 
significantly neutral, this may be a hint of large cosmic 
variance in the epoch of reionization. Indeed, a strongly 
fluctuating t^s is itself a signature of the pre-overlap 
era. In the extreme interpretation that the Lya and Ly/3 
troughs of SDSS J1030+0524 indicate xm ~ 0.2 down 
to z — 5.95, this implies large modulation in the ioniza- 
tion fraction and typical bubble sizes of order Az = 0.38, 
or L ~ 150 Mpc comoving. This is substantially larger 
than the H II regions of these extremely bright and rare 
quasars, i?Hii ^ 30 Mpc comoving. It is also larger than 
theoretical expectations for the scale of typical H II re- 
gions at the tail end of reionization (Furlanetto et al. 
2004; Wyithe & Loeb 2004b), so we consider such a sce- 
nario to be extremely unlikely. At the other extreme, 
the measured Toff are compatible with a slightly faster 
increase in Toff along the sightlinc to SDSS J1030+0524 
(Toff > 9.9, 2cr), than to SDSS J1148-t-5251 , (Toff « 6- 
10). Whether a highly-ionized universe can tolerate such 
scatter is unclear. We have shown that a simple Toff 
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analysis is a blunt instrument, and more sophisticated 
interpretation is required to put strong constraints on 
reionization. 

Note that surveys of Lya emitters have found no evolu- 
tion of the luminosity function of Lya emitters at z = 5.7 
and z = 6.5 (Malhotra & Rhoads 2004, Stern et al 2004), 
which has been interpreted as evidence against percola- 
tion taking place at z 6, in agreement with our con- 
clusions. 

What other signatures of large Lya optical depth could 
emerge from the spectra? One possibility is the O I ab- 
sorption forest (Oh 2002). O I is an excellent tracer of 
neutral hydrogen: it has a very similar ionization po- 
tential E = 13.62eV, so it lies in tight charge exchange 
equilibrium. Furthermore, its 1302 A absorption line lies 
redward of the Lya forest. Although the O I forest coin- 
cides with a noisy portion of the night sky, the sightline 
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